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Abstract 

Solutions of the Dirac equation with spin and pseudospin symmetry for the scalar 
and vector trigonometric scarf potential in D-dimensions within the framework of an 
approximation scheme to the centrifugal barrier are obtained. The energy spectral 
and the two-component spinor eigenfunctions are obtained. The bound state energy 
levels for various values of dimension D, n, k and the potential range parameter a 
are also presented. 
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1 Introduction 

In other to investigate nuclear shell structure, the study of spin and pseudospin symmetric 
solutions of the Dirac equation has been an important area of research in nuclear physics. 
The concept of spin and pseudospin symmetry with nuclear shell model has been used 
widely in explaining a number of phenomena in nuclear physics and related area. 

Within the framework of the Dirac equation, the spin symmetry occurs when the dif- 
ference of the potential between the repulsive Lorentz vector potential V(r) and attractive 
Lorentz scalar potential S(r) is a constant, that is, A(r) = V{r) — S{r) = constant. 
However, pseudospin symmetry arises when sum of the potential of the repulsive Lorentz 
vector potential V(r) and attractive Lorentz scalar potential S(r) is a constant, that is, 
E(r) = V{r) + S{r) = constant [H El S H El il [7] . 
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The spin symmetry is relevant for meson [8j . The pseudospin symmetry has been used: 
to study the structure of the deformed nuclei [9l [10], to construct an effective shell model 
coupling scheme [11] and identical bands and triaxality observed in nuclei pj]- 

Recently many researchers have applied spin and pseudospin symmetry conditions on 
numbers of potentials. These potentials include: Hulthen potential [131 [E] , Eckart poten- 
tial [m [ISl [13 US], Poschl-teller potential [H [2D], the Rosen- Morse potential ^ [22], 
three parameter diatomic potential [IH], harmonic potential [23], Pseudoharmonic poten- 
tial [21], Manning- Rosen potential [25], Wood-Saxon potential [26] and kratzer potential 
with angle dependent potential [27] . 

On the other hand, exact analytical solutions for i states are not possible for some 
of these potentials due to the centrifugal term. An approximation to the (pseudo or) - 
centrifugal term has been used by many researchers to obtain the approximate analytical 
solutions [m [151 [ISl [2H1 [221 [23- Several researchers have used various methods to solve 
spin and pseudopsin symmetry problems including the centrifugal approximation rang- 
ing from the Asymptotic Iteration method (AIM), the Nikiforov-Uvarov method (N-U), 
Supersymmetric and shape invariance method to functional analysis method. Due to its 
wider applications and simplicity , we are adopting the Nikiforov-Uvarov method [301 E] • 
This method find its usefullness in the study of relativistic and non-relativistic quantum 
mechanics, see Ikhdair et al. (2011) and the references therein [32] . 

In this study, we are considering the Scarf potential [33] 

V{r) = (1) 
sm ar 

where Vq is the depth of the Scarf potential and a is the range of the potential. 

Futhermore, with interest in higher dimensional spaces, the multidimensional relativis- 
tic and non-relativistic equations have been investigated by many authors. This investi- 
gation have employed Isotropic harmonic oscillator plus inverse quadratic potential [31], 
Pseudoharmonic potential [3S] , Kratzer- Fues potential [3S1 [3Z] and hydrogen atom [3S] (for 
comprehensive reviews on abtrary dimension, see Dong 2007). In the present work. Section 
2 contained a brief review of the Nikiforov-Uvarov (N-U) method. The Dirac equation in 
D-dimension is described in Section 3. Section 4 contains the bound state solutions of 
the Dirac equation with the trigonometric Scarf potential and conclusion is contained in 



Section 5. 



2 The Nikiforov-Uvarov (N-U) Method 

The N-U method is based on solving a second order hnear differential equation by reducing 
it to a generalized equation of hypergeometric type [301 EI]- By introducing an appropriate 
coordinate transformation z = z{r), this equation can be re- written in the following form. 



aiz) 



a'^iz) 



(2) 



where (7{z) and a{z) are polynomials, at most of second degree, and t{z) is a first degree 
polynomial. Now, if one takes the following factorization 



^{z) = y{z)(t){z). 



The equation ([3]) reduces to a hypergeometric type equation 



(3) 



a{z)y"{z)+T{z)y\z) + Xy{z) = Q, 



(4) 



where 

t{z) =f{z) + 2Ti{z). (5) 
The function tt and the parameter A require for this method are define as follows 



a'[z] -Tlz 



TTlZ 



± 



'a'{z) 



TiZ 



a{z) + ka{z) 



and 



A = k + 7i'(z). 



(6) 



(7) 



In other to find the value of k, the expression under the square root must be a square 
of a polynomial. This gives the polynomial nlz) which is dependent on the transformation 
function z{r). Also the parameter A defined in equation ([7]) takes the form 



A = A. 



— UrT — 



Uririr - 1) „ 

cr 



(8) 



The polynomial solutions yn{z) are given by the Rodrigue relation 



VnAZ) 



p{z) dz^ 
3 



(9) 



where A^^^ is the normahzation constant and p{z) is the weight function satisfying 

^Ja{z)p{z)]=T{z)p{z). (10) 
Second part of the wave function (j){z) can be obtained from the relation 

7r(.) = a(.)^[ln</>(z)]. (11) 

3 The Dirac Equation in D- Dimensions 

The D-dimension Dirac equation for a particle of mass M, with radially symmetric Lorentz 
vector V{r) and Lorentz scalar potential S{r) is given(in atomic units h = c = 1 [201 HOI 

i—^ = H^{f), H = J2 (^jPj +/3[M + S{r)] + V{r), (12) 
at .^^ 

where M is the mass of the particle, V{r) and S{r) are spherically symmetric vector and 

scalar potentials, respectively, P is the momentum operator, {ttj} and /3 are Dirac matrices 

satisfying anti-commutation relations. 

Following the procedure stated by Ciftci et al. (2005) and Agboola (2011), the spinor 

wave functions can be written using the Pauli-Dirac representation as: 

where F„,,K(r) and Gn,.K{r) are the radial wave functions of the upper- and the lower - 
spinor components, respectively. Y-^iyLo) and Y-^iVLr)) are the hyperspherical harmonic 
functions coupled with the total angular momentum, while, £ and I correspond to the 
orbital and pseudo-orbital angular momentum quantum numbers for spin and pseudospin 
symmetry conditions, n = — {j + "^^) fo^' aligned spin j = (■ + \ and n = {j + "^^) fo^" 
unaligned spin j = (- — \. 

Substituting Equation f|T3|) into Equation f|T2|) . the two coupled second-order differential 
equation for the upper and lower component spinors are obtained follows: 

+ -)j FnUr) = [M + + S{r) - V{r)] G„.,.(r), (14) 
+ 7) Gn^Ar) = [M- + S{r) + V{r)] F„„,«(r) . (15) 



T 2 



/ D-2\ 
«: = ±(j + ^— ) (13) 



Substituting the expression of ^^^^^(r) obtained from equation f|T5|) into equation f|T4|) . we 
have two second-order differential equations for the upper and lower component spinors 
as: 



+ [M + - A(r)] [M - + S(r)] F„,,.(r) 



dA{r) / d _|_ K\ 
dr ^dr r ' 

[M + En^,, - A(r)] 



(16) 



k,(k — 1] 



dr'" 



+ 



+ [M + ^„^,« - A(r)] [M - + S(r)] G„,,.(r) 



dr ydr r ' 



(17) 



where A(r) = y(r) — S{r) and S(r) = y(r) + S{r) are the difference and the sum of the 
potentials V{r) and 5(r), respectively and k = ± ( ^^+^"^ ) [201 SHI US] • 

In the presence of the spin symmetry condition, that is, the difference potential A(r) = 
V{r) — S{r) = Cs = constant or = 0, then, equation f|T6|) reduces into 



/«(/«+ 11 



+ 



+ [M + E„^,«-C,]S(r) F„,,«(r) 



£;^^,,-M^ + C.(M-i?„,,.) 



Then, the lower component Gnr^l^) Dirac spinor is obtained as 

1 



d ^ K 
dr r 



M + - a 

where En^^^^ + M 7^ 0, only real positive energy state exist when = 

Again, under the pseudospin symmetry condition, that is, the sum potential S(r 
V{r) + S{r) = Cps constant or = 0, then, equation (fT7|l can be written as 



(19) 



d"^ k{k — I] 



dr' 



+ 



[M-K.,. + Cp,]A(r) G„,,4r) 



_M2 + Cp,(M-K,,,) G„,,,(r 



(20) 



where the upper component Fn^^i^{r) is obtained as 

1 



FnJr) 



Mc^ - E^^ + 
5 



d 
dr 



(21) 



where Enn — Mc ^ 0, only real negative energy state exist when Cpg = [22| 

These equations cannot be solved exactly due to centrifugal term, so, we used the 
following approximation 

1 

- ^ (22) 

r"^ cos"^ ar 

With this approximation, the second order differential equations in equation (ITS]) and (120]) 
can be solved approximately by using Nikiforov-Uvarov method. 

4 Bound State Solutions of the Trigonometric Scarf 
potential 

4.1 Spin symmetry solutions of the Dirac equation with the 
Trigonometry Scarf potential with arbitrary k in Y> dimen- 
sion 

Under the condition of the spin symmetry, A(r) = V{r) — S{r) = Cg, the vector and scalar 
trigonometric Scarf potential is considered [33]. Then, 

V(r) = ^ and S{r) = 1^ (23) 

sin ar sin ar 

where Vo and 5*0 are the depth of vector and scalar trigonometric Scarf potential, therefore, 

EW = (24) 
sm ar 

Substituting equations ( 122]) and (124]) into equation (|T8]) . we have 

72 „,^,2 ^^,2 



7" s^a^ + ^ 



sm ar 



Fn^,k{r) = 0, (25) 



where 



-f = k{k + 1), 

£2^2 = - E„,.,,) (M + - Cs) , 

rya2 ^ ^ ^ ^^^^^ _ ^ (26) 



In other to apply N-U method, we introduced a new variable 

z = — tan^ ar 
6 



(27) 



With equation f l27|) . equation f l25|) becomes 
1 - 32 „, _ 1 



-7^' + (7 + - r?)^ + . (28) 



nr,k\~; • 2^(1 -z)^"-"^-^ ' 4^2(1-2)2 

Substituting for a{z),a{z) and t{z) in equation ([6]), we obtained function ni^z) as 



'1 + 477) + AieJ{l + 477) - 4£2 



2 + 877 + 4z£j(l + 4r/) 



+ (1+477). 

(29) 

According to the Nikiforov-Uvarov method, the expression under the square root must be 
a square of a polynomial and therefore, the new 7r(2;) function for each k is obtained as 



[A + {2ie - 1)] + 1 [1 - A] ; A; = -i [// - 7 - _ 
-f [A + (2ze-l)] + i[l + A];fc = -i[r/-7-e^]- 



Me 



Aie 



Aie 



[A - {2ie + 1)] + i [1 - A] ; A; = -i [77 - 7 - 
-f [A-(2ze-l)] + l[l + A];A; = -i[7/-7-£2]+-^ 



(30) 



where A = y/1 + Ar]. 

For the polynomial t{z) = f{z) + 27r , with a negative derivative, we select 

n{z) = -|[A + {2te - 1)] + ^[1 + A]; k = -^Ir^ - j - e'] - ^ 

and when this is combined with A = A; + tt'Iz), gives 

A = 

Now using equation ([8]), we obtained 



(31) 



(32) 



A = —271^. + 2ienr + Auj. + 2nr{nr — 1), 
and on comparing equations fl33|) and fl32|) . we obtained the energy equation as 



(33) 



£ = - [2(271, + 1) + (2A; + 1) +A] 



(34) 



Or more explicitly 

(M-^„,,,)(M + K.,.-C, 



a 
'7 



2 r 



.n 2 



271 + D + - + (4v; + So){M + - 



where we have defined a principal quantum number as ti — £ = 2ri.r + 1, chosen k = 



(35) 



and 7z = 1, 2, .... The explicit values of the energy for different values of 7^,, a and D 
are shown for spin symmetry in the Table 1. below: 
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Table 1: Bound-state energy level for different values of n^, a and D with M = Cg = 



D 


rij. 


a = 0.0001 


a = 0.0005 


a = 0.001 


a = 0.005 


a = 0.01 




1 
1 














2 


3.0006 


3.0029 


3.0057 


3.0285 


3.0570 


O 




^ 0007 


3 0037 


3 0073 


3 0367 


3 0733 




4 


3 OOOQ 


3 0045 


3 OOQO 


3 0448 


3 08Q5 




(J 


3 001 1 


3 0053 


3 0106 


3 0350 


3 1057 




1 


o.uuuo 


^ 0094 


^ 004Q 


^ 094 


^ 04SQ 




2 


3.0007 


3.0033 


3.0065 


3.0326 


3.0652 


4 


3 


3.0008 


3.0041 


3.0082 


3.0408 


3.0814 




4 


3.0010 


3.0049 


3.0092 


3.0489 


3.0976 




5 


3.0011 


3.0057 


3.0114 


3.0570 


3.1138 




1 


3.0006 


3.0029 


3.0057 


3.0285 


3.0570 




2 


3.0007 


3.0037 


3.0073 


3.0367 


3.0733 


5 


3 


3.0009 


3.0045 


3.0090 


3.0448 


3.0895 




4 


3.0011 


3.0053 


3.0106 


3.0530 


3.1057 




5 


3.0012 


3.0061 


3.0122 


3.0611 


3.1219 




1 


3.0007 


3.0033 


3.0065 


3.0326 


3.0652 




2 


3.0008 


3.0041 


3.0082 


3.0408 


3.0814 


6 


3 


3.0010 


3.0049 


3.0098 


3.0489 


3.0976 




4 


3.0011 


3.0057 


3.0114 


3.0570 


3.1138 




5 


3.0013 


3.0065 


3.0131 


3.0651 


3.1300 
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From the Nikiforov-Uvarov method, the hypergeometric function y{z) depend on the 
determination of the weight function p{z) which is satisfied by the differential equation 



dz 



[a{z)p{z)] = t{z)p{z). Thus p{z) is obtained as 

p{z) = {i-zy^z^, 



by substituting p{z) into the Rodrigue relation (j9]) one gets 



Urir 



Z) ^Z 2 



T'^il-z 



(36) 



(37) 



The solution of equation (1371) may be expressed in terms of the Jacobi polynomials of the 
form [S] 

yn. = Nn^pt^\l-2z) (38) 

and (t>(z) is obtained as 



l+A 



z) = 1 - zY'z^ 



(39) 



therefore, the radial wavefunction of the upper component of the Dirac wavefunction 

Fnr,k{z) can be written as 



Fn^A^) = iVn.(l - zfz'^Pt'\l - 2Z). 



(40) 



The corresponding lower -component spinor Gnr,k{z) in equation f|T9l) can be obtained 
by using the following relation [15] : 



z . 



(41) 



Using above equation pTl) . we obtained the lower component spinor Gnr,k{z) from equation 
(HH]) as 



Gnr,k{z) 



iV„„(l-^)f.^p£^'^^(l-2z) 



+ iV„^(l-2;) — P„V-i nl-2^) 



•.A + l\/l-z\ ^ k 
-OL < I z — - — — ^ + eJz - 



z tanh a 



A + 2nr 
-a <, I I z I + e 



where Nn,. is the normalization factor to be determined from normalization condition 

(43) 



/ I F„^,fc(r) |2 dr = -— / Fn^4z)[^{l - z)fdz = 1. 
Jo lia Jo 



This can further be written as 



-2ia Jo 



z^(l — z 



iie— 1 



1 2 



dz. 



(44) 



Two different forms of tlie Jacobi polynomials are ^5] HU 1^ : 
and 



p=0 



" n\r{n + c + d+l)^,\r 



T{n + c + d + r + 1) fz~l 
r + c + 1 



(45) 



(46) 



where 



n 



r(n+l) 



J- y r!(n-r)! r(r+l)r(n-r+l) 
_ . ., , ,2 

We obtained 



Pjf'^^^l - 2z) from equations and (jlH]) as 
P^'^)(l - 2z)f = z^^'^-^il - zYAr.,{p, r) 
and putting equation (H7I) into equation (H^ . we have 

—2ia Jo 

Using the definition of hypergeometric function 

z'^-^l - z)-^ds = -[ 2Fi(a, b;a + l; 1)], 
Jo a 

we obtained N„^ as 



(47) 



(4J 



where 



-l)"[r(ri + c+l)]2r(n + rf+l) 



(49) 



(50) 



E 



'-l)P+'q''~P+'T{n + c + d + r + l) 



T{n + c- p + l)T{n + d+ l)T{c + d+ 1) fr^Q p\r\{n - p)\{n - r)\T{p + d + 1) ' 



Bnq{p, r) 



2Fi(i[2 + A] + n ~ p + r,l - is - p; ^[4: + A]n - p + r; 1) 



-ia{A +[2n-p + r + l]) 



(51) 



we have used c = ie, and d = ^. 



4.2 Pseudopin symmetry solutions of the Dirac equation with 
the Trigonometry Scarf potential with arbitrary k: in D di- 



mension 



In the case of pseudospin symmetry, that is, the sum I](^) = V{r) + S{r) = Cps, we have 

(K - So) 



E( 



sm ar 



(52) 
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Substituting equations fl22|) and f l52|) into equation fj20|) . the following equation is ob- 
tained, 

i^n.,fc(r) = 0, (53) 



2 2 

TpQ^ 2 2 _|_ Vp'^ 



(ir^ cos^ ar ' sm" ar 
where, we have used the following notations: 

7p = k{k- 1), 
52^2 = (M + E„„,,) (M - + Cps) , 

r^ptt^ = - So) (M - E^^^k + C,) . 



(54) 



Using the same procedure (the N-U method)used in solving equation (^5^, the energy 
equation is obtained as follows: 



(M+^„,,«)(M-^„,,,+C, 



ps J 



a 



2n + D-2 + -^a2 + (4K - 5o)(M - + C^,) 



(55) 

where we have defined a principal quantum number as n — i = 2nr + 1, chosen k = 1+ '■'^"^^ 
and n = 1,2, .... The explicit values of the energy for different values of n^, a and D 
are shown for pseudospin symmetry in the Table 2. below: 

We can consequently obtain the wave function ^ and G^r./t' respectively as 



(56) 



1 -2z 



GnrAz) 



+ a^X^ - z)^ z— Pl_, 



i tanh a 



-1,^+1 



A„ + 2n, 
-a <; 1 1 — - ) + Ep 



(57) 



where Ap = y/1 + rjp. 

The normalization constant C„^ is obtained as 



where 
Dngip,r) 



^ Enqip, r) 
-l)"[r(n + c+ l)]2r(n + (i+ 1) 



(5^ 



(-l)p+'-g"-p+T(n + c + (i + r + l) 
r(n + c-p+ l)r(n + (i+ l)r(c + d+ 1) ^ j9!r!(n-p)!(n-r)!r(p + (i+l) ' 



Enqip, r) 



'2Fi{l[2 + Ap] + n-p + r,l-iep-p;l[A + Ap]n-p + r;l] 



1 -1 



-za(Ap+ [2n-p + r + l]) 



(59) 



we have used c = and 



2 • 
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Table 2: Bound-state energy level for different values of n,., a and D with M 



D 


rij. 


a = 0.0001 


a = 0.0005 


a = 0.001 


a = 0.005 


a = 0.01 




1 
1 




-z.uuuuuuoo 




9 nnnn^^^^ 

-^.yjijyjyjoooo 


9 nnni ^^^9 

-Z.UUU-LOOOZ 




2 


-2.00000003 


-2.00000075 


-2.00000300 


-2.00007500 


-2.00030000 


O 




-9 ononono^ 


-9 000001 


-9 OOOOO^S'^ 


-9 0001 S^S'^ 


-9 0005^394 




4 


-9 onononos 


-9 00000908 


-9 OOOOOSS'^ 


-9 000908^9 


-9 00083310 




(J 


-9 onononi 9 


-9 00000^00 


-9 00001900 


-9 000^0000 


-9 001 1QQ^9 




1 


9 00000009 


-9 0000001^9 


-9 0000090?^ 


-9 0000'=i90R 


-9 00090839 




2 


-2.00000004 


-2.00000100 


-2.00000408 


-2.00010207 


-2.00040828 


4 


3 


-2.00000006 


-2.00000168 


-2.00000675 


-2.00016874 


-2.00067484 




4 


-2.00000010 


-2.00000252 


-2.00000101 


-2.00025206 


-2.00100799 




5 


-2.00000014 


-2.00000352 


-2.00000141 


-2.00035204 


-2.00140767 




1 


-2.00000003 


-2.00000075 


-2.00000300 


-2.00007499 


-2.00030000 




2 


-2.00000005 


-2.00000133 


-2.00000533 


-2.00013333 


-2.00053324 


5 


3 


-2.00000008 


-2.00000208 


-2.00000833 


-2.00020832 


-2.00083310 




4 


-2.00000012 


-2.00000252 


-2.00001200 


-2.00003000 


-2.00119952 




5 


-2.00000016 


-2.00000408 


-2.00001633 


-2.00040828 


-2.00163245 




1 


-2.00000004 


-2.00000102 


-2.00000408 


-2.00010208 


-2.00040828 




2 


-2.00000007 


-2.00000168 


-2.00000674 


-2.00016874 


-2.00067485 


6 


3 


-2.00000010 


-2.00000252 


-2.00001008 


-2.00025206 


-2.00100799 




4 


-2.00000014 


-2.00000352 


-2.00001408 


-2.00035204 


-2.00140767 




6 


-2.00000019 


-2.00000469 


-2.00001874 


-2.00046868 


-2.00187383 
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4.3 Special cases 

4.3.1 PT Symmetric trigonometric Scarf potential 

We consider a case in which a of equation ([1]) as a pure imaginary parameter, that is 
a — )■ ia, then the potential becomes [33] 

So-Vo 



V{r) 



sinh ar 



(60) 



For the spin and pseudospin symmetry conditions, the corresponding energy eigenvalues 
equation are obtained, respectively as. 



a 



(M - ^„„.)(M + -Cs) = — 



2n + D- -jA{So + Vo) (M + - C,) - 
a * 



(61) 



and 



(M+E„,,«)(M-K,,«+Q 



a 



ps J 



2n + D-2- -jA{So - K) (M - K.,. + Cp,) - 



(62) 



4.3.2 PT Symmetric and g-deformed trigonometric Scarf potential 



In equation ([T]), we have used a mapping r — )• r — Mnyg[33l HT], then, the potential 
becomes 



VJr 



sinhg ar 



(63) 



and the energy equations for spin and pseudospin symmetry of the g-deformed trigono- 
metric Scarf potential are obtained, respectively as: 



(M -^„^,,)(M + -C) 



and 



{M+En^,n){M-En^^n + Cps) 



a 



2 r 



2n + D- -JAq-^So + K)(M + - C,) - a' 
a * 



a 
T 



(64) 



2r2 + D - 2 - -^4g-i(.So - K)(M - K.,. + Cps) - 



(65) 



4.3.3 Non-PT Symmetric, non Hermitian and g-deformed trigonometric Scarf 
potential 

Another form of the potential in equation ([T]) is obtained by putting Sq = Si + iS2 and 
Vo = Vi + iV2, a = ia and q ^ iq [33]. In this case, after some manipulations, the potential 
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takes the form 



V,„=4(^4^)±^(^. (66) 

[e^ — iqe 



Then, the energy equations of the spin and pseudospin symmetry for this potential are 
obtained, respectively as: 



{M-En^,,){M+Er,^,,-Cs) 

2n + D- 1^4g-i[(yi + Si)i + (S2 + V2)](M + - C,) - 



(67) 



and 



a 



2n + D-2- -^Aq-^[(Vi - Si)i + (S2 - V2)](M - E^^,, + C^s) - 

5 Conclusion 



2 



(68) 



In this work, we have obtained the bound state solutions of the £)-dimensional Dirac 
equation with spin and pseudospin symmetry for scalar and vector trigonometric Scarf 
potential. The two -component spinors and corresponding the energy equation have been 
obtained using Nikiforov-Uvarov method. The numerical results shows that the bound 
state energy level E^^^,^ for spin and pseudospin symmetry increases with both dimensions 
D and range parameter a. In addition, PT symmetric, g-deformed, non-PT symmetric 
and non-Hermitian version of this potential are investigated in this study as special cases 
and their energy equations are obtained. 
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